clc;
syms r x y; 
f=x*y/(x^2+y^2); 
L=limit(subs(f,y,r*x),x,0) 